Two dimensional electron liquid in the presence of Rashba spin-orbit coupling: 
symmetric momentum space occupation states 
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The orientation of the local electron spin quantization axis in momentum space is identified as 
the most significant physical variable in determining the states of a two-dimensional electron liquid 
in the presence of Rashba spin-orbit coupling. Within mean field theory several phases can be 
identified that are characterized by a simple symmetric momentum space occupation. The problem 
admits uniform paramagnetic as well as spin polarized chiral solutions. The latter have a nontrivial 
spin texture in momentum space and are constructed out of states that are not solutions of the non 
interacting Hamiltonian. The concept of generalized chirality as well as the stability of spatially 
homogeneous states are also discussed. 
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PACS numbers: 71.10.-w, 71.70.Ej, 71.45. Gm, 75.30.Fv 

Recent interest in the physics of two dimensional elec- 
tronic systems subject to the effects of Rashba spin-orbit 
interaction makes the corresponding many-body problem 
one of timely interest. Since it is well known that inter- 
action effects in modern quasi-two dimensional electron 
liquid devices are quite strong in a wide range of accessi- 
ble densitieS)i»2i2i£ it is clearly of fundamental importance 
to establish a sound set of theoretical notions about the 
effects of the electron-electron interactions in this prob- 
lem. Patterning the general approach on that followed 
for the familiar case of the electron liquid in the absence 
of spin-orbit^ the first step in tackling the many-electron 
problem is to establish a meaningful mean field theory. 
Once this step is successfully taken, and the relevant dy- 
namical variables and symmetries identified, correlation 
corrections can be eventually approximately calculated 
by a variety of methods. The purpose of this paper is to 
describe some of the general results of the mean field the- 
ory of the two dimensional electron liquid in the presence 
of Rashba spin-orbit^ The present analysis is limited to 
the simple case of symmetric momentum space occupa- 
tion states that, as we will show, can be completely and 
elegantly classified in terms of what we will refer to as 
the generalized chirality. 

The problem is defined by the following model hamil- 
tonian: 



» = E^ + ^E 



(1) 



where the electronic motion is limited to the x — y plane 
and the single particle terms contain a spin-orbit inter- 
action of the Rashba type££ 



H = l^ + a (ZxPy - OyPx) 



(2) 



where a is assumed positive and terms corresponding 
to an homogeneous neutralizing background are under- 
stood. In the homogeneous case Hq has eigenfunctions 



and eigenvalues given by 



Vk,±(r) = 



±1 

j p i<f>k 



h 2 k 2 



=F ah k 



2m 

(3) 

where L is the linear size of the system and tf>k is the angle 
between the direction of the wave vector and the x-axis. 
Interestingly spin-orbit forces each state in momentum 
space to have its own spin quantization axis. This di- 
rection lies in the x — y plane and makes an angle of ? 
with k. The corresponding unit vector will be denoted 
by ^k- These states form two split bands characterized 
by opposite chirality which are schematized in Fig. [T] 7 

Since the mean field theory involves, but as we shall see 
is not limited to, single Slater determinants obtained by 
occupying the single particle states (jHJl, it proves neces- 
sary to establish a complete yet manageable classification 
scheme. We have found that for states with isotropic, 
compact occupation in momentum space the following 
quantity, 




for < xo < 1 
for Xo = 1 



(4) 



which wc will refer to as the generalized chirality, offers 



such possibility. In |gj) xo = 



JV+-JV- 



is the ordinary 



N + +N- 

chirality as defined in terms of the occupation N± of each 
band, while ki n , k out and k± are geometrical parameters 
in momentum space characterizing the occupied regions 
(circles for the isotropic case) and are defined in Fig. [T] 
Depending on the relative occupation of the bands, x can 
acquire any value from zero to infinity and allows one to 
uniquely label the relevant set of many-particle states for 
fixed density n = N/L 2 . This must be contrasted with 
the ordinary chirality which satisfies the condition xo < 1 
and equals unity for all the states in which only one of 
the chiral bands is occupied. 

Limiting for the time being our analysis to spatially 
homogeneous solutions, it is readily established that, be- 
cause of the peculiar physics of the Rashba spin-orbit, 
the momentum space local orientation of the electron 
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FIG. 1: Two distinct, symmetric compact ways of occupying 
non interacting chiral states in k space. In the top panel 
only the lowest chiral band is occupied and the generalized 
chirality differs from the ordinary chirality which in this case 
is unity. The left panel represents the non interacting energy 
spectrum. 




FIG. 2: Geometry of the local (in momentum space) spin 
quantization axis. The z-axis is taken perpendicular to the 
plane of motion. 



spin quantization axis, as determined by its unit vec- 
tor Sk, plays a crucial physical role. In Fig. [2] the angles 
/3k and 7k determining the orientation of Sk are defined. 
Accordingly our idea is to construct many-particle Fock 
states "J [rik ± , Sk] with occupations nk ± using single par- 
ticle plane wave states each characterized by its own Sk- 
Our main task is therefore reduced to determine what, 
for a given density, is the most energetically favorable set 
of nk± and corresponding spin quantization axes Sk- 

From a formal point of view this problem can be 
handled by making use of a standard Wick decoupling 
procedure^ We find however more instructive and ele- 
gant to minimize the total energy as a functional of both 
rik± and Sk- It is a simple exercise to write down the 
expectation value of the full hamiltonian of the electron 
liquid ([!]) taken over a generic single Slater determinant 
*[nk±,Sk]: 



E 

k; fi=± 



2u2 



2m 



- E Uk - k ' (- 1 + w' ^ k ' * k ') "kc^V -( 5 ) 



E 

k,k'; [i,/j,'=i: 



The corresponding, renormalizcd, single particle energies 
£k^ can then be readily obtained by differentiation with 
respect to n\ SL ±^ 

We begin our analysis of Eq. ([5]) by noting that, as 
one can show, for symmetric occupations, the case 7k = 
0k + § corresponds to an energy minimum. In this case 
by implementing the momentum rescaling 



/3k = /3(- 



0, 



(6) 



the functional can be expressed in Rydberg units as 
follows 



[nk,u, Sk] 



a- 



(7) 



where we have also introduced the dimensionless param- 
eters rj 1 — 'it «bti and a = This formula is quite 
remarkable for it displays a simple explicit dependence 
on the spin-orbit coupling constant and an high degree 
of universality. Specifically IC is a simple universal func- 
tion of x- 



IC(x) = 0(1 - X )(l + X 2 ) + 0(x-l)2 X 



(8) 



while 1Z and £ x are in addition universal functionals of 
the rescaled momentum dependent angle [3 whose expres- 
sions are reported elsewhere^ For paramagnetic states 
the situation simplifies further since for this class of solu- 
tions /3k = P — § and only the dependence on x survives. 
Figure [3] displays the typical dependence of 1Z and £ x on 
X for actual solutions /3 of the problem for different val- 
ues of a. As we show below, for a given value of the 
latter, the final values of 1Z and £ x can be expressed as a 
function of x on ly- Notice the cuspy behavior about the 
critical value \ — 1 which separates states for which one 
or two (renormalized) chiral bands are occupied. 

At this point the mean field energy can be obtained by 
minimization with respect to both x an d Pk- 



£MF(a,r s 



(9) 



For a fixed x minimization with respect to /3k reduces 
to finding the lowest lying solutions of the corresponding 
mean field equation 



tan /3(re) 



=d9 + Aita k 



2 — 2kk' cos t 



■ d6 



here simplified to the relevant case 7k 
sponding to a vanishing in plane polarization. Remark- 
ably this equation is devoid of any explicit dependence 
on r s , the solution /3 being only dependent on a and x- 
Notice that (3 = 5 is always a solution of Eq. (jTUJ) . The 
latter has rather interesting symmetries. For instance 
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FIG. 3: Plot of K\x, P] (top) and £ x [x, P] (bottom), as denned 
in ((7J), evaluated for the mean field solutions determined by 
solving Eq. (|10[) as functions of the generalized chirality \ f° r 
different values of a. Dashed lines: unpolarized case. Solid 
lines: transverse ferromagnetic solutions. 
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FIG. 4: Top: universal spin azimuthal angle /3(k) for the case 
of a = 0.1 and various values of \- Bottom: same for the 
corresponding reciprocal values of x- The dots mark the end- 
ing points of the [yjl — xj , \A + Xj intervals corresponding 
to the occupied regions of Fig. [T] The top (bottom) panel 
corresponds to the top (bottom) panel of Fig. Q] 




one can see that /3i/ x (k) = ^(y^K). Eq. (fT0|) makes it 
manifest how, for this class of solutions, the generalized 
chirality x is the central parameter of the problem. As an 
illustration, typical solutions for the universal azimuthal 
angle (3(n) are displayed in Fig. [4] There, solutions with 
j3 ^ 77 correspond to polarized states in which the spin 
polarization points along the z-axis. Interestingly, these 
are constructed out of single particle states that are not 
solutions of the single particle hamiltonian ([2]) and dis- 
play an intriguing spin texture in momentum space. This 
is shown in Fig. [5] where three possible cases are depicted. 

Further minimization with respect to \ allows one to 
compare the energies of these homogeneous, symmetri- 
cally occupied states. This leads to the determination of 
the boundaries (lines of first order transitions) between 
the paramagnetic (PM) and the z-axis polarized phases 
(FZ, dotted area) in the relative phase diagram of Fig. [5] 
Within the ferromagnetic region the generalized chiral- 
ity is constant and equals one. This phenomenon is due 
to the cuspidal behavior of £ x for this value of \- The 
fractional polarization is readily seen to be given by 

2 (S z ) rVt+x 
p{a,x) = t. at = K cos/3(«)dK . (11) 

p vanishes for unpolarized states (J3 = 5) and acquires 
an r s independent value in the polarized phase for which 




FIG. 5: Momentum space spin texture for polarized states. 
Top: x — 0-4 and x = 2.5 for a — 0.1. Bottom: Lowest 
energy polarized state for a — 0.2 (with x — 1 and maximal 
polarization) . 



P !J. In general the ferromagnetic polarization is less 
than one as it can be surmised from Fig. [7] On the other 
hand p — ► 1 for vanishing a as one expects to recover 
the fully polarized Bloch ferromagnetic state of the two 
dimensional electron liquid in the absence of spin-orbit 
coupling. 5 A most remarkable feature of this phase dia- 
gram is the reentrant paramagnetic phase at lower den- 
sities, something due to the density dependence of the 
Rashba term (see Eq. 0). 

We discuss next the possibility of inhomogeneous 
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FIG. 6: Relative phase diagram in the {r s ,a) plane of the 
isotropically occupied, spatially homogeneous mean field so- 
lutions. In the shaded region the gas is polarized along the 
perpendicular direction. Along the dashed line the energy of 
the unpolarized states is minimized by x = !• 




FIG. 7: Fractional polarization for the transverse ferromag- 
netic solutions (with x = 1) 5 plotted as function of a. 



where we have chosen Q = 2kp£ = —x as to con- 
nect states on opposite sides of the Fermi energy and 
suitable for the case of €f = in Fig. Q] The corre- 
sponding leading change in the total mean field energy 
is a functional A£j./f [^] of the linear variations of 
the matrix elements of the single particle density matrix 



dla/3 



$) whose expression can be found in 



Pa/3 , .. , 

Ref. [5|. At this point, after some algebra, one finds that 
the following judicious choice of the amplitudes 



(bk F )2 |rt k±Q -7t k | 



I sin 6 



|k±3|i 



| cos <t> k± Q I 







bk F < |k ± S| < e bk F 
otherwise, 



( 13 ) 

where b is an arbitrarily small positive number, leads to 
a negative AEmf even without allowing for momentum 
space repopulation. 

An inspection of our procedure reveals that, within 
mean field theory, several of the effects of the interac- 
tions can be described in terms of the concept of effective 
magnetic field, a quantity defined as 



9VB 



B, 



-ha k i 



2L 2 ^ 

k' 



(nk'+-nk'-)«k-k' Sk' 



(14) 

the last term being associated with the exchange energy. 
As it can be shown, this concept can be employed to 
write the mean field equation (|10|) at once. 



states. It is clear that in the low density limit the sys- 
tem will form a Wigner crystal. On the other hand, we 
have found that, as in the absence of spin-orbit coupling, 
within the mean field approximation, homogeneous states 
are unstable to spin density wave type distortions^ We 
have been able to construct a rigorous proof of this state- 
ment valid for all densities. The procedure is patterned 
after that of Ref. [f| One assumes the following set of 
trial single particle wave functions: 

V>k(r) ~ y>k,+ (r) + A k+ <p k+Q . + (r) + A k _</? k _ Q .+ (r) , 

(12) 



We conclude by noticing that, as in the absence of spin- 
orbit, one should reasonably expect correlation effects to 
change the density range in which the various transitions 
occur. On the other hand the qualitative structure of the 
phase diagram should remain essentially unchanged. We 
also expect that the interplay of interactions and spin- 
orbit will lead to interesting novel phenomena not only 
in the transition region to the Wigner crystal, where the 
spin polarization of the states will play a major role, but 
also at intermediate densities where interesting in-plane 
spin structures could be stabilized by both exchange and 
correlations. The current analysis only offers a glimpse 
at the rich physical complexity of this problem. 
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